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Abstract 

We reexamine relic abundance of a singlet fermion as a CDM candidate, which contributes 
to the neutrino mass generation through radiative seesaw mechanism. We search solutions 
for Yukawa couplings and the mass spectrum of relevant fields to explain neutrino oscil- 
lation data. For such solutions, we show that an abundance of a lightest singlet fermion 
can be consistent with WMAP data without conflicting with both bounds of fi — * cy 
and r — > fi'y. This reconciliation does not need any modification of the original radiative 
seesaw model other than by specifying flavor structure of Yukawa couplings and taking 
account of coannihilation effects. 
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1 Introduction 



Observation of nonzero neutrino masses [T] and the existence of dark matter [2] gives 
big impact on the study of physics beyond the standard model (SM). Elucidation of the 
origin of neutrino masses and dark matter is now one of the biggest subjects in this 
field. Evidence for small neutrino masses suggests some additional structure in the SM 
as ones required in the seesaw mechanism [3j. On the other hand, in order to explain 
the abundance of dark matter, we need some symmetry to guarantee the stability of 
a dark matter candidate. The most famous example of such symmetry is R parity in 
supersymmetric models. Even in non-supersymmetric models in general, if Z 2 symmetry 
is required by some natural reason, it can play the same role as R parity as long as the 
ordinary SM fields have its even charge. From this point of view, there appears a very 
interesting idea that neutrino masses may be intimately related to the existence of cold 
dark matter. 

If we consider that neutrino masses are generated radiatively under the assumption 
that a model has Z 2 symmetry whose existence is justified to forbid tree- level Dirac neu- 
trino masses pE], this symmetry can guarantee the stability of a Z 2 odd neutral particle 
which may be cold dark matter (CDM). In this direction the relation between neutrino 
masses and CDM has been studied in various articles [5j EJ [8] . Since such scenarios re- 
quire the introduction of new particles and interactions in the SM, however, it can induce 
dangerous effects in various phenomena as usual. In fact, it is suggested that contradic- 
tion could appear between the strength of Yukawa couplings to satisfy the required relic 
abundance of CDM and the bound for fi — > e'j as long as a singlet fermion is considered 
as a CDM candidate [6]. Some attempts have been proposed to overcome this fault by 
modifying the model [9]. 

In this paper we fix our target on a minimal model in [4], which can generate neutrino 
mass radiatively. And we also confine our study into the case that a CDM candidate is 
one of the singlet fermions. We reanalyze whether its CDM abundance can be consistent 
with lepton flavor violating processes only by specifying detailed structure of the neutrino 
mass matrix and taking account of coannihilation effects. Our result will show that the 
CDM abundance can be consistent with lepton flavor violating processes even within the 
simplest radiative seesaw framework. 

The paper is organized as follows. In the next section we define our model briefly 
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and discuss neutrino masses and mixing for the explanation of neutrino oscillation data 
based on this model. In section 3 we show our result on both relic abundance of the CDM 
candidate and constraints from the lepton flavor violating processes. We summarize the 
paper in section 4. 

2 Neutrino mass due to radiative effects 

We consider a model which is an extension of the SM with an additional SU(2)x doublet 
scalar rj and three gauge singlet right-handed fermions N k @]. The model is also imposed 
by Z 2 symmetry to forbid tree-level Dirac masses for neutrinos. We assign odd charge of 
this Z 2 symmetry to all of these new fields, although Z 2 even charge is assigned to all of 
the SM contents. 

Lagrangian relevant to N k invariant under the SM gauge symmetry and this Z 2 sym- 
metry are written as 

C N = (mrfdpP R N k ) + ~ {M k W k P R N k + M;W k P L N c k ) - (hak&iP R N k + h.c), (1) 

where £ a stands for a lepton doublet and a charged lepton mass matrix is assumed to 
be diagonalized. We note that N k can have mass terms invariant under the imposed 
symmetry. For simplicity, these masses M k and Yukawa couplings h ak are assumed to be 
real in the following discussion. 

Scalar doublets $ and rj have invariant scalar potential 

Zi Zi 

+ ^A 5 [(^r/) 2 + h.c], (2) 

where $ is the ordinary SM Higgs doublet. If we assume that only $ obtains a vacuum 

expectation value (VEV) such as ($°) = v but rj obtains no VEV, neutrinos cannot have 

tree-level Dirac masses. However, neutrino masses can be generated radiatively through 

a one-loop diagram which has rj° and N k in internal lines. This radiative masses can be 

small as long as A5 is sufficiently smalll^ Since A5 is assumed to be very small, masses of 

real and imaginary parts of rj and also ^ are considered to be degenerate and they can 

§ Since we can introduce a new U(l) symmetry in case of A5 = 0, the smallness of A5 may be considered 
as a natural assumption. 
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be written as m 2 , = m 2 + (A3 + A4)?; 2 . Thus, in the following discussion we use this Wq as 
the mass of 77. 

Radiatively generated neutrino masses are expressed by using the Yukawa couplings 
h a k and three mass scales A k as 

3 

(•M-v)a/3 — 2^ hakhpk-h-ki (3) 



k=i 



where A k is defined by 



A. 



87T 2 M fc V mo 



Kx) 



x 



1-x 2 



1 + 



1 -x 2 



In a;' 



(4) 



By using this mass matrix, we now consider how to explain neutrino oscillation data. 
Since it is known that neutrino oscillation data are well explained by using the Maki- 
Nakagawa- Sakata (MNS) matrix^ 



U 



( cos 6 sin 6 \ 

sin 8 cos 6 1 

V2 V2 y/2 

sin 9 cos 1 

\ y/2 V2 V2 / 



(5) 



we assume that the mass matrix (Ej) is diagonalized as \J T M. V XJ = diag(mi, m,2, m^) 
Then, we find that the following diagonalization conditions should be satisfied: 
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£ (2/4 sin 26 + 2^2h ek {h^ - h rk ) cos 26 - (h rk - h^) 2 sin 2fl) = 0, 

k=l 

3 3 

h e k (h^k + h rk ) = 0, — /i T fc) (h^ + h rk 



0. 



(6) 



k=l k=l 

The mass eigenvalues are expressed as 

3 

nil 



^ ^efe cos2 61 + ^= sin ZOKkiKk - h^ k ) + i sin 2 (/i Tfe - h^ k f^j A k , 

m 2 = ^ (hlk sin2 ® ~ sin 26h ek (h rk - h^k) + ^ cos 2 # (/i rfe - /i Mfc ) 2 ^ A k , 
3 1 

m 3 = ^2 9 ^ rfc + h^k) 2 A k . 



(7) 



fc=i 



^In this matrix sin# 13 is assumed to be zero, for simplicity. Since sin# 13 is expected to be very small, 
the present analysis is considered to be straightforwardly extended to the case with sinf^ ^ 0. 
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Here, among various solutions for the conditions (EJ), we consider a simple solution such 

as 



h f 



0, h^i h T i, h e j ^ 0, hfj,j 



(9) 



where i ^ j is assumed ^\ This means either of i or j runs two values of k — 1, 2, 3 such as 
i — 1,2 and j = 3, for example. By substituting this in the first condition in ([6]) we have 

+ a 1 he i 
tan 6/ = = — -. 

The mass eigenvalues (j7|) are rewritten as 

m i = f^ej cos # + V2h T j sin 6> j A 3 - = 

m 2 = (h e j sin — \[2h T j cos 6 1 j Aj 



2h 2 - 

l Jhi a 

cos 2 ' 



m 3 = 2^.A i5 (10) 

where the summation for i and j should be understood. We use eq. (jHD in the last equality 
for m li2 t 

Now we impose phenomenological requirements on the model. If we recall that the 
MNS matrix is given by U defined in eq. (J5]), it is found that we can use sin 2 (9 ~ 0.33 
and m 2 ~ 7.66 x 10~ 5 eV 2 suggested by the solar neutrino and KamLAND data and also 
777.3 — 2-46 x 10~ 3 eV 2 suggested by the atmospheric neutrino and K2K data PQ. As this 



result, we obtain 



h 2 Tj Aj ~ 2.9 x 10~ 3 eV, 



hl-Ai ~ 2.5 x 10" 2 eV, 



'ID 



where the summation on i and j is abbreviated. Similar one-loop diagrams to the one for 
neutrino masses contribute to the lepton flavor violating processes like £ a — > Ibl- It gives 
the most severe constraint on the model. Its branching ratio is estimated as [TT] 

3 

Br(£ a -> 47) 



y2hi ak h e k F 2 ( — - J 



where Fz(x) is given by 



Fo(x) 



647r(G F m 2 ) 2 
1 - 6x 2 + 3x 4 + 2x 6 - 6x 4 In x 2 



(12) 



6(1 -x 



2\4 



(13) 



" If = —h T k is satisfied for all k, it is also a solution for the diagonalization conditions. However, 

such a solution cannot satisfy neutrino oscillation data and then we do not consider this case here. 
**This type of neutrino mass hierarchy induced from the mass matrix ([3]) has been considered to analyze 

neutrino oscillation data in other context [TO] , 
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The present upper bounds for Br([i — > e7) and Br{r — > fij) are given as 1.2 x 10 11 [12] 
and 6.8 x 10~ 8 [13], respectively. If we use eq. (jSj), these constraints can be written as 



If we discuss other phenomenological features of the model, we should analyze them under 
the conditions ffTTl) and (JHJ). 

To proceed with the study of CDM abundance in the next section, it is convenient to 
classify possible cases for the relation between Mk and m . Since we assume N\ to be a 
CDM candidate and we are interested in the effect of coannihilation, N\ is considered to 
be almost degenerate with other Z 2 odd fields. In such a situation, physically distinctive 
cases may be classified as 

(i) M 1 < M 2 < M 3 , mo, (ii) M x < M 2 , m < M 3 , (iii) M x < m < M 2 , M 3 . 

Although there is no logical correlation between Yukawa couplings and masses of singlet 
fermions without introducing some symmetry, here we only assume that the masses of the 
corresponding singlet fermions are equal if Yukawa couplings are equal in eq. (jSJ) • We can 
identify important processes for determination of CDM abundance under this assumption. 
If we take % = 1, 2 and j = 3, two possible cases (i) and (ii) should be considered. In case 
(i) we need to take account of coannihilation between Ni and N 2 , in which only Yukawa 
couplings are relevant to this process. On the other hand, in case (ii) we take account of 
coannihilation among N\, N 2 and rj. Gauge interaction is expected to play an important 
role in this case. If we take i — 1 and j = 2, 3, the case (iii) with M 2 = M3 is a target 
for the investigation. In this case coannihilation between N\ and 77 is expected to play a 
crucial role. Both gauge and Yukawa interactions are relevant to this case. Although final 
states of coannihilation in the cases (ii) and (iii) can include antiproton, the case (i) can 
not include it but include only lepton pairs. This aspect makes the case (i) interesting in 
the relation to the PAMELA e+ and p data [H], and also the ATIC/PPB-BETS (e+ + e~) 
data [151 IB] • We will come back to this point later. In the next section we will confine 
our study to this case. Other cases will be discussed elsewhere. 




(14) 
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3 Coannihilation of the CDM candidate 



In this section we consider (co) annihilation of N\ through Yukawa couplings in the case (i). 
For the estimation of the relic abundance of Ni, we follow the method given in [17], which 
is developed to take account of coannihilation effects. If we introduce the dimensionless 
parameter x as x = Mi/T, the decoupling temperature Tf of N% can be estimated by 
using effective cross section a eS and effective degrees of freedom g cS as 

Xf - m 1/2 1/2 ' ^ i£) ' 

9* x f 

where v re \ is the relative velocity of annihilating fields. a e g and g e s are defined as 

2 2 

<7eff = —a NlNl + 2 — 2 — °N 1 N 2 { i - + A) e +— %aJV 2 (l + ^j e > 
9 c s 9 c s 9 c s 



9cs 



g Nl +g N2 (l + Af 2 e- Ax , (16) 



where m pl = 1.22 x 10 19 GeV and A is defined by A = (M 2 - M 1 )/M 1 . If we define a eff 
and b e ff by a e g \v Te \\ = a e fr + b e «v^ el , thermally averaged cross section can be written as 
(c" e ff |f re i|) = a c s + 66 e fj/x. In the following analysis, A ~ is assumed since we consider 
the case (i). Thus, if we use this decoupling temperature Xf, the relic abundance can be 
estimated by 

!»• = m im X ^ . (17) 

gj m p i(GeV)(a eff + 3b eS /x f ) 

(Co)annihilation proceeds via t-channel exchange of rf and 77 through Yukawa inter- 
actions. The final states are composed of only leptons fi, r, u^, v T and their antiparticles. 
We note that antiproton is never produced. The (co)annililation cross section of and 



N i2 is estimated as 



Ml 

~ 2\2 



1 Ml{M\- 
12tt (Aff + mg) 



3(Mf + mlf rcl 



1 M?(M* + mg) 2 v-, 

1 0-TT- f )U2_|_ rr ,2'S4 ^rel 2 ; Qil Q ' 2 ^ ^ 2 ' 



a,/3 

where ii,Z2 should be considered as 1 or 2. As can be seen from this expression with 
H = H — 1) the annihilation of iVj occurs only through a p-wave channel. On the 
other hand, coannihilation defined by i\ ^ 12 can have s-wave contributions in general. 
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Fig. 1 Regions sandwiched by red solid lines satisfy the WMAP data flh 2 = 0.11 ± 0.06 for CDM 
abundance. Blue dashed and blue dash-dotted lines show the bounds for /i — ► and r — > fjtj, respectively. 
Green dotted lines represent contours for 6a ^ — l.Ox 10~ n . The values of r\ and r% are fixed as r\ = 0.8 
in both graphs and r% — 4 and 10 in the left and right graphs. 

However, if we take account of the assumed conditions jSJ) in our model, we find that s- 
wave contributions cancel out and only p-wave contributions remain. Thus, the relevant 
cross section is found to be written as 

1 Ai?(M« + mg) 2 
^ ilNi2 M = - {M!+ml y KiAAv (19) 



We have a e fr = and 

,2 , ,2 \2 



127im 2 {1 + rlY 

where we define r\ by r\ = Mi/m . Applying this effective cross section to eq. (|T7|) . we 
can estimate the relic abundance of Aq. 

In the case (i), the constraints (11 ip obtained from the neutrino oscillation data are 
written as 

h 2 r3 A 3 ~ 2.9 x KT 3 eV, (h 2 rl + h 2 r2 )A 1 ~ 2.5 x KT 2 eV. (21) 



If we use this relation in eq. ( 1201) . the effective annihilation cross section of Aq can be 
expressed by using A x instead of Yukawa couplings h T \ and h r2 . As found from eq. (THl) . the 
bound of — > e7 directly constrains only /i t3 which is not relevant to the N% annihilation 
in this case. Since the relevant Yukawa couplings h T ± and h T 2 are constrained by the 
bound of r — > /ry, this bound can be much more severe constraint than that of /i — > ej. 
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mo 



mo 



Fig. 2 The same figures as Fig. 1 for a different r\. Although n is fixed as t\ = 0.4 in both graphs, rs 
is put as 4 and 10 in the left and right graphs as in Fig. 1. 



These constraints can also be written as conditions on Ai as follows, 



Ai > 3.0 



r 3 /(ri) 



nl(r 3 ) 

A 1 > ( 0.34F 2 (n 



F 2 (r 3 ) 



500 GeV\ : 



m 



eV, 



4.0 x 10- 2 ^iHF 2 (r 3 ; 



500 GeV^ 



m 



eV, 



(22) 



ri/(r 3 ) 

where we define r 3 by r 3 = M 3 /mo- We use a relation A 3 = (ri/(r 3 )/r 3 /(ri))Ai obtained 
from eq. (j3J) and also the conditions in eq. (|2T|) in this derivation. 

Using these results, we plot favorable regions in the (mo, Ai) plane by fixing the values 
of ri and r 3 . Since our considering case (i) corresponds to r% < 1 and r 3 > 1, we fix these 
to some typical values. We plot examples of the allowed regions for (n, r 3 ) = (0.8, 4) and 
(0.8, 10) in Fig. 1 and also for (ri, r 3 ) = (0.4, 4) and (0.4, 10) in Fig. 2. In both figures, blue 
dashed and blue dash-dotted lines show the bounds for // — > and r — > fij, respectively. 
The upper regions of both lines satisfy these constraints. Thin bands sandwiched by the 
red solid lines corresponds to region to realize the value of flh 2 required by the WMAP 
data. Form this figure we find that the relic abundance of N x can be consistent with the 
WMAP data without conflicting the bounds of lepton flavor violating processes. Figs. 
1 and 2 show that mass of N\ should be larger than 700 GeV and 500 GeV for each 
case. A3/A1 takes values of 0(1), for example, 1.14 and 0.80 for r 3 = 4 and 10 in case 
of r% = 0.8, respectively. Although A x has rather small values ~ 0.01 in the allowed 
regions, Yukawa couplings h T \ and h r2 can be confirmed to be in the perturbative regions 
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by taking account of eq. ( I2T1) . Since these couplings contribute tor-> /i7, this constraint 
can be much stronger than fi — > e'j as found from Figs. 1 and 2. However, we can find 
consistent solutions by fixing r% and r 3 suitably. 

The reconciliation between the CDM abundance and the lepton flavor violating neutral 
processes can be shown to be accomplished even in the original radiative seesaw model 
without substantial modification of the model. In the present flavor structure of Yukawa 
couplings, the (co) annihilation of N\ and fi — > ej are induced by the different ones, 
respectively. Yukawa couplings relevant to the (co) annihilation of Ni contribute to r — > 7/7 
whose bound is much weaker than \i — > e^y. This feature makes their reconciliation 
possible by arranging the masses of singlet fermions so as to satisfy the requirements from 
the neutrino oscillation data. 

Finally we give remarks on some predictions of the model. Both direct and indirect 
detections of dark matter are crucial to judge whether the considering model for dark 
matter is viable or not [THf [T9] . Since the CDM candidate has couplings only to [i and r, it 
can decay to these. Model independent analysis of the data of PAMELA and ATIC/PPB- 
BETS has been done in [20l EI]- Ref.[2T] suggests that the best fit is obtained for M ~ 
1 TeV with CDM annihilating into and a good fit is obtained for M ~ 2 TeV with 

CDM annihilating into t + t~. It is interesting that this is consistent with our results 
obtained in the present analysis. Encouraged by this result, we would like to add some 
qualitative arguments on the related subjects in our particular model. 

As shown in eq. (1191) . N\ annihilation cross section a\v\ is dominated by p-wave con- 
tribution due to helicity suppression. However, although p-wave contribution which has 
v 2 dependence dominates the annihilation cross section at freeze out time where v ~ 0.2, 
it is largely suppressed in the present Galaxy where v ~ 10~ 3 . This makes s-wave contri- 
bution relevant to PAMELA anomaly rather than the p-wave contribution. Since s-wave 
annihilation cross section can be estimated as cr\v\ ~ it 1 wttt z\2 where rxit is the mass 
of final fermions, N\ annihilation in the Galaxy occurs mainly through N\Ni — > t + t~. 
If we use typical values of m , h T \ and r 1 obtained as the solutions consistent with the 
WMAP data in this paper, we find that the boost factor should be O(10 6 ) or larger HU The 
model cannot induce this amount of enhancement for the annihilation cross section in the 
present form. For the explanation of this boost factor, there may be two possibilities: (i) 



ft We note that a\v\ ~ 10 23 cm 3 /sec is required to explain the positron excess in the PAMELA data. 
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the model should be extended such that the relic Ni has a large non-thermal component as 
discussed in |22j, for example, or (ii) it should be explained by some astrophysical effects. 
However, it seems difficult to obtain substantial effects by a simple extension referred in 
case (i). Since we can not make the annihilation cross section of N% itself larger preserving 
the features of the model, main effect should come from the increase of number density of 
the non-thermal component of relic N% as a result of the decay of other fields. However, 
it is severely constrained by the WMAP data and we have no freedom to obtain the large 
boost factor mentioned above. 

Bremsstrahlung from the charged fields associated to this annihilation N±Ni — > t + t~ 
yields diffuse photons. We may check the model by comparing the flux of diffuse photon 
expected from this N\ annihilation based on the PAMELA data with observations such as 
Hess and Fermi/GLAST. Such a model independent analysis is presented in [23]. Although 
it suggests that dark matter annihilation with r + r _ final states may be difficult to be 
consistent with diffuse photon data, the assumptions in that analysis seems not to be 
applied to our model. On the other hand, dark matter annihilation into final states 
composed of three fields such as e + e~7 can be dominant processes if the annihilation cross 
section is helicity suppressed [21] . The radio emission from synchrotron radiation and 7- 
ray emission from inverse Compton scattering from the charged fields produced by the 
dark matter annihilation may be also useful to discriminate the origin of positron excess 
|25j . Observational data of diffuse photon obtained in the Fermi/GLAST experiment may 
give us crucial hints for these [26]. Anyway, detailed analysis of diffuse photon is necessary 
to check the validity of the present model. 

Decay of r ± also produces neutrino flux. If we use the PAMELA positron data, its 
flux can be roughly estimated as O(10 5 ~ 6 ) GeV/(cm 2 -sec-str) at the relevant neutrino 
energy. This flux is larger than neutrino flux expected from certain types of AGN but 
smaller than the atmospheric neutrino flux [2Z]- This suggests that the N% annihilation 
is difficult to be detected through the observation of the neutrino flux on the Earth. 

Direct detection of Ni is also an interesting subject. Although 77 has the interac- 
tion shown in eq. (pQ) and no direct interaction with quarks, Ni can be scattered by nuclei 
through one-loop effect with Z boson exchange. Since Ni is a Majorana fermion, this effec- 
tive interaction with quark is expressed by an axial vector interaction d q N\ 757^ A"i 5757^5 

2^2 ^,2 

with d q ~ u n )2^i T 3q , which yields spin dependent scattering. If we use the parame- 
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ters obtained in this paper, this spin dependent cross section is roughly estimated as 
0(1CT 41 ) cm 2 . This is much smaller than the present bound of spin dependent elastic 
scattering cross section for dark matter with 0(1) TeV mass [28]. Thus, it seems difficult 
to find this dark matter even in the next generation direct detection experiments. Studies 
related to these aspects of the similar model can also be found in [29] although lepton 
number violating constraints are not taken into account there. 

In addition to these indirect and direct search of dark matter, there may be some other 
phenomena which could show characteristic features of the model. The effective mass in 
the neutrinoless double j3 decay is given as a fixed value m e g = Am s 2 ol sin 6 so \ ~ 2.9 x 1CT 3 
eV, which is one order of magnitude below the reach of near future experiments. Present 
data for the magnetic dipole moment of muon shows discrepancy between a value predicted 
by the SM and the experimental result |30j. In the present model there is one-loop 
contribution to 5a M , which is estimated as [TT] 



This shows that these values predicted by our model are two orders of magnitude smaller 
than Sa^ = (30.2 ± 8.7) x 10~ 10 [31J. In order to improve this situation for the 5a^, 
additional contributions to 5a^ are required in our model. Such contributions may be 
obtained by embedding our scenario in supersymmetric models, in which an ordinary 
supersymmetric CDM candidate such as the lightest neutralino does not dominate the 
required relic abundance. We will discuss such extensions elsewhere. 

4 Summary 

The radiative seesaw model considered in this paper is one of interesting possibilities to 
explain the origin of neutrino masses. It can include a cold dark matter candidate as 
an important ingredient of the neutrino mass generation. However, the model has been 
considered to have a severe discrepancy between magnitude of Yukawa couplings required 
by the dark matter relic abundance and the suppression of lepton flavor violating neutral 




(23) 



A contour for 5a^ 



1.0 x 10 is also plotted by a green dotted line in the figures. 
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processes. In this study we have proposed a new possibility to relax this tension within 
the original minimal radiative seesaw model without introducing additional interactions. 
We have found that the model can overcome this problem simultaneously satisfying the 
conditions required by the neutrino oscillation data as long as the Yukawa couplings and 
also the mass hierarchy of the singlet fermions have appropriate structure. The present 
study shows that even the minimal radiative seesaw model can be an interesting candidate 
for models which relates neutrino masses to the existence of dark matter. 

The model may be relevant to the PAMELA e + and p data, and also the ATIC/PPB- 
BETS (e + + e~) data. We have briefly presented qualitative observations on the detection 
of the diffuse particles produced in the dark matter annihilation and the direct search 
of this dark matter. Since these analyses are rough and qualitative ones, we need more 
quantitative study to mention on the predictions of the model in detail. This scenario 
may play an important role in some supersymmetric models if it is embedded in the 
supersymmetric framework. Although we have considered only the model with restricted 
coannihilation processes here, other cases are also expected to give interesting possibilities. 
These points may be worthy for further study and will be discussed elsewhere. 

This work is partially supported by a Grant-in-Aid for Scientific Research (C) from Japan 
Society for Promotion of Science (No. 17540246). 
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